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We investigate tlie dynamical spin and cliarge structure factors and the one-particle spectral function of the 
one-dimensional extended Hubbard model at half band-filling using the dynamical density-matrix renormaliza- 
tion group method. The influence of the model parameters on these frequency- and momentum-resolved dynam- 
ical correlation functions is discussed in detail for the Mott-insulating regime. We find quantitative agreement 
between our numerical results and experiments for the optical conductivity, resonant inelastic X-ray scattering, 
neutron scattering, and angle-resolved photoemission spectroscopy in the quasi-one-dimensional Mott insulator 
SrCu02. 
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I. INTRODUCTION 

The spectral properties of low-dimensional correlated elec- 
tron systems have attracted a significant amount of interest 
in recent years.' In particular, the energy and momentum 
resolution of spectroscopic techniques have been greatly im- 
proved. These methods now provide plenty of information 
about the dispersion and intensity of electronic excitations in 
low-dimensional materials. For instance, the dynamical sepa- 
ration of charge and spin degrees of freedom, which is generic 
in the theory of quasi-one-dimensional correlated systems, 
has been observed directly in angle resolved photoemission 
spectroscopy (ARPES) of strongly anisotropic cuprate com- 
pounds;-'' From the viewpoint of theory, however, calculating 
the response functions that are probed in these experiments 
remains a great challenge and thus a coherent quantitative de- 
scription of these materials is still lacking. 

Few exact analytical results are available for dynamical re- 
sponse functions in itinerant correlated electron systems. In 
gapless systems dynamical correlation functions can be eval- 
uated analytically with field-theoretical methods in the limit 
of small excitation energiesil Gapped systems can also be 
treated analytically in the weak-coupling limit when the gap 
is small compared to the electron band-widthi^ At finite ex- 
citation energies exact analytical results are available in the 
strong-coupling regimai^ or when additional symmetries are 
presentsSi 

There are several numerical approaches for calculating 
momentum- and frequency-dependent correlation functions 
in lattice models of coiTelated electrons. These include Ex- 
act DiagonalizationpS Cluster Perturbation Theory (CPT)|ii 
Variational CPTfi^ and Quantum Monte Carlo algorithms li^iil 
Here we employ the Dynamical Density-Matrix Renormal- 
ization Group (DDMRG) method.Ji It is a zero-temperature 
method which has been successfully applied to the study of 
dynamical properties in metallic and insulating phases of one- 
dimensional extended Hubbard model s>i^iiii^ii2i2£ 

In this paper we consider the dynamical properties of the 
half-filled extended Hubbard model with nearest-neighbor in- 
teraction. It is believed to be a minimal model for quasi-one- 



dimensional Mott-insulators, such as the cuprate compounds 
SrCu02 and Sr2Cu03. In a previous worlii^ we have shown 
that this model describes the optical conductivity and the dis- 
persive structures observed in the resonant inelastic x-ray scat- 
tering (RIXS) experiments of SrCu02. Here we show that the 
model also captures essential results obtained in inelastic neu- 
tron scatteringri electron energy loss spectroscopy (EELS)p^ 
and ARPES^"^ experiments. To this end we present numeri- 
cal results for the energy- and momentum-resolved spin and 
charge structure factors and the one-particle spectral function 
of the extended Hubbard model. We also discuss in detail the 
influence of the model parameters on the line shapes of struc- 
ture factors and spectral functions. 



II. MODEL AND METHODS 

The one-dimensional extended Hubbard model is defined 
by the Hamiltonian 



H 



-UYl \ "'-T -oil 



(1) 



- 1 



- 1 



Here c\ ^, c-i ^ are fermionic creation and annihilation oper- 
ators for a particle with spin a =t, i at site / = 1, . . . ,L, 

In the following we 
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consider chains with open boundary conditions and an even 
number of lattice sites. The total number of electrons is equal 
to the number of lattice sites (half-filled band). 

For V/t = the model reduces to the simple Hubbard 
Hamiltonian. Without magnetic field the physical excitations 
of the insulating Hubbard model are combinations of collec- 
tive modes, spinons and (anti-)holons, for any U > 0^ The 
spinons are gapless and charge neutral excitations that carry 
spin 1/2. Anti-holons and holons are gapped modes with 
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charge ±e and no spin. When the nearest-neighbor interac- 
tion V is turned on and U > 0, the system remains a Mott- 
insulator until V ^ U/2. For larger V the system becomes a 
charge-density wave (CDW) insulator>ai In this work we only 
consider the Mott-insulating regime U > 2V > 0. The pa- 
rameters U = 7.8t, V = 1.3t and t = 0.435eV^ were shown in 
Ref. 02 to describe both the optical conductivity and the dis- 
persive structures observed in the RIXS spectrum for SrCu02. 
In addition these parameters yield an effective spin-exchange 
coupling J Ri 0.24cV which is in good agreement with recent 
high-precision neutron scattering data.Si 

The dynamical spin structure factors S{q,uj) is the imagi- 
nary part of the spin-spin correlation function 

S{q,Lj) = -lm{^o\S~^^^ ^— 5,1^o), (2) 

where jV'o) and Eq are the ground state wave function and 
energy and 77 ^ 0+. The operator is the Fourier transform 
of the local z-component of the spin operator S[ = ni,^ — ni,i. 
We set ?i = 1 everywhere, so that an excitation momentum 
and energy are equal to the wavevector q and the frequency 
Lo, respectively. When we replace by the Fourier transform 
fiq of the local particle density ni = hi^^ + fii^i we obtain the 
dynamical charge structure factor 

N{q,Lj) ^ -lm{ijo\nl-^ 5-- figl^o). (3) 

TT H + LU - Eo-trj 

Similarly, the one-particle spectral function is defined through 

A{q,uj) = -lm{^Jo\cl,^ ^— -c (4) 

71" H + uj - Ea - IT] 

where ^ is the Fourier transform of the annihilation opera- 
tor C;^. 

These dynamical correlation functions can be calculated 
with the DDMRG method for a finite system size L and a 
finite broadening 7].-^ Therefore, all DDMRG spectra pre- 
sented here are a convolution of the true dynamical correlation 
function with a Lorentzian distribution of width rj. The spec- 
tral properties in the thermodynamic limit can be determined 
using a finite-size scaling analysis with an appropriate size- 
dependent broadening 77 (L) cx 1/L. In general, momentum 
dependent operators = S^,hq, or ^ are defined from lo- 
cal operators a; = , , or ^ using the one-electron eigen- 
states of the non-interacting system with periodic boundary 
conditions, i.e. 

aq = l=Y,e-''i'ai (5) 

with momentum q = 27rZ/ L (we set the lattice constant equal 
to 1) and integers -L/2 < Z < L/2. Since the DMRG algo- 
rithm performs best for open boundary conditions, we extend 
the definitions of S{q, lo), N{q, uj), and A{q, uj) and write 



aq = J X! sin(A:Oa/ (6) 
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FIG. 1: Line shapes of the spin structure factor S{q,uj) for < 
q < TV calculated with a broadening rj — 0.05t for U — 7.8t and 

V = 1.3t. 

with (quasi-)momenta q = tiZ/{L + 1) for integers \ < Z < 
L. The expansion coefficients are the eigenstates of a free par- 
ticle in a box. Both definitions of aq should become equivalent 
in the thermodynamic limit L ^ 00. Comparisons with Bethe 
Ansatz results for the Hubbard model have confirmed that 
this procedure gives accurate results for g-dependent quan- 
tjtJesJiiM2, In our DDMRG calculations we have used up 
to 400 density matrix eigenstates. For all results presented 
here truncation errors are negligible compared to the finite 
resolution in momentum [Aq = tt/{L + 1)] and in energy 
{Alli ^ 7] (X l/L) imposed by finite system sizes. 

III. DYNAMICAL SPIN STRUCTURE FACTOR 

In this section we present our results for the dynamical spin 
structure factor S{q,u!). We show that the spin dynamics 
of the extended Hubbard model agrees with a recent high- 
precision inelastic neutron scattering experiment in SrCu02i^ 
In addition, we discuss the influence of the model parameters 
on the spin structure factor of the Hamiltonian 

Figure[2shows the line shapes of the dynamical spin struc- 
ture factor in the extended Hubbard model for U ~ 7.8t and 

V = 1.3<. The system size is L = 100 and we have used 
a broadening 77 = 5t/L. The spectrum is dominated by the 
spectral weight at the lower onset which disperses from zero 
at q = 0, TT to a maximum at 7r/2. Above the low energy peak 
there is a continuum of spectral weight which is bounded by 
an upper onset. 

These line shapes are strongly reminiscent of the spin struc- 
ture factor of the spin- 1/2 Heisenberg chain, which in good ap- 
proximation is described by the Miiller-AnsataSi. The Miiller- 
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FIG. 2: (a) Density plot of the data in Fig.Q (b) Upper (circles) and 
lower (squares) onset of the two-spinon continuum extracted from 
the DDMRG data. Solid lines are fits to the dCP dispersion relations. 
The dotted line is the dCP dispersion using the exchange coupling 
Jexp = 0.226 eV from neutron scattering experiments in SrCu02!~ 

Ansatz structure factor is given by 

Q (n,A , e{uj-uji^{q))e{uju{q)~uj) 
^MA(g,wj = A , „ „, , = (/) 

where A is a prefactor and wu.lI?) are the des Cloiseaux- 
Pearson (dCP) dispersion relations. They describe the com- 
pact support of the two-spinon continuum of the spin- 1/2 
Heisenberg model 

ttJ 

^L{q) = — I sin(g)| , (8a) 
ujviq) = ^J|sin(g/2)| . (8b) 

where J is the nearest-neighbor spin-exchange coupling. The 
Miiller-Ansatz structure factor is non-zero only within the 
bounds of ujuxiq) and there is a square-root divergence at 
the low-energy onset. 

The similarity of SuAiq,'^) with S{q,uj) in the extended 
Hubbard model can be seen more clearly in a density plot. 
Figure 13 a) shows a density plot of the spectral weight distri- 
bution. Most spectral weight is located at the lower onset of 
the two-spinon continuum as in the Heisenberg model. Fig- 
ure |2jb) shows the lower and upper boundaries of the two- 
spinon continuum extracted from the DDMRG data. The up- 
per boundary was obtained by analyzing the second derivative 
d'^S{q,u!)/dLu'^ and the lower onset was identified with the 
position of the low-energy peak. These boundaries are very 
similar to the dCP-dispersion relations dSal and j8bt of the 
Heisenberg model. Fitting our numerical data to these rela- 
tions we obtain effective exchange couplings J ~ 0.242cV 



FIG. 3: Dependence of the spin structure factor at momentum g ~ tt 
on U and V. The broadening is = Ws/ 20 in the left panel and 
ri = in the other two panels. 



(lower boundary) and J = 0.248eV (upper boundary), re- 
spectively (using t = 0.435 eV). These values are in good 
agreement with the value J = 0.24eV reported in our pre- 
vious worbi^ and with the coupling Jexp = 0.226eV deter- 
mined from inelastic neutron scattering data for SrCu02i^ 
As the Miiller-Ansatz Q with the exchange coupling Jexp = 
0.226cV describes the neutron scattering spectrum of SrCu02 
very well|2i we conclude that the experimental spectrum is 
also well described by the spin structure factor of the extended 
Hubbard model for U = 7.8t, V = 1.3t and t = 0.435 eV. 

A recent study^ of the spin structure factor in the Hubbard 
model has shown that there are significant itinerancy effects 
[a transfer of spectral weight in S{q,u!) due to the coupling 
of the spin excitations to charge fluctuations] at low and inter- 
mediate values of the Hubbard interaction U. In this regime 
the spin structure factor of the Hubbard model differs from the 
spectrum 5'ma(<7, t^) of the Heisenberg model (i.e., the strong- 
coupling limit U ^ t of the Hubbard model), where charge 
fluctuations are suppressed and electrons are completely lo- 
calized. To illustrate this effect we show DDMRG spectra 
obtained for various parameters U and V on 60-site lattices 
in Fig. |3l We have normalized S{q ~ 7r,Lj) by the to- 
tal spectral weight ^(Tr) = J S{Tr,uj)duj, which we obtain 
by calculating the ground state expectation value S{7t) = 
(V'o|S'q=,r'5'q=:^|V'o)- In thcsc figurcs the energies are given 
in units of the two-spinon band width 2Ws, which is equal to 
IT J in the Heisenberg model. For V/t = the value of Ws is 
known exactly from the thermodynamic Bethe-Ansatz solu- 
tionis^ When V/t we determine Wg from the upper onset 
of the dynamical spin structure factor 

Figure |3j a) shows S{7t,uj) calculated with a broadening 
r/ = 1^^/20 for V/t = and three values of U/t = 0, 1, 3. 
At U/t = we recover the case of free electrons and ac- 
cordingly most spectral weight is found in the peak at the up- 
per onset {uj = 2Ws) of S{q, uj). When U/t becomes larger 
the spectral weight is quickly redistributed towards the lower 
onset at = 0. Already at U/t = 3 the upper peak is no 
longer visible and most spectral weight is found at the lower 
onset. This supports the analysis of Ref. but shows that 
itinerancy effects become rapidly weaker for increasing U. 



4 



a) U/t=7.8, V/t=0 b) U/t=7.8, V/t=3 
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FIG. 4: Spin structure factor ^(g, lv) at high energies for (a) U/t = 
7.8, V/t = and (b) U/t = 7.8, V/t = 3. The broadening is 
rj = 0.2t. Vertical lines indicate the one-particle gaps Ac. 
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FIG. 5: Density plot of the renormalized charge structure factor 
N{q,u)/q^ for U/t = 7.8, V/t = 1.3. For g > 7r/2 a solid line 
follows the spinon dispersion (shifted by the value of the charge gap) 
and triangles indicate the position of the spectrum onset in the nu- 
merical data. 



We now set V/t = 1.3 and vary U/t. The resulting curves 
S{Tr,Lu) are shown in Figure |3jb). At U/t ~ 3 the largest 
part of the spectral weight is at the low-energy onset. With 
increasing U /t the low-energy peak becomes more and more 
pronounced until only little spectral weight remains at the up- 
per onset. For larger interaction strengths the change of the 
next-nearest neighbor interaction has a weaker influence on 
the spectral distribution. Figure|3jc) shows 5'(7r,cj) for con- 
stant U/t = 7.8 and varying V/t. While the weight of the 
peak at = is clearly diminished with growing V, the ef- 
fect appears to be less pronounced than in the case of smaller 
U . In summary, our DDMRG results for the extended Hub- 
bard model confirms that itinerancy effects are important in 
the dynamical spin structure factor of Mott insulatorsi^ We 
note, however, that despite the coupling between spin and 
charge sectors the upper and lower onsets of the two-spinon 
continuum are very well described by the dCP relations 
and dSal of the Heisenberg model from strong to intermediate 
interaction as illustrated by the results in Fig.|2| 

Our numerically results reveals another effect of charge 
fluctuations. The total spectral weight of the spin structure 
factor is not exhausted by low-energy excitations ijj{q) < 
uj\j{q). There are dispersive features in that spectrum at en- 
ergies above the one-particle gap, Ac = £'o(l) — i?o(l) + 
2i?o(0), where Eo{x) denotes the ground state energy of the 
Hamiltonian Q with x electrons added to or removed from 
the half-filled system. Figure |3 shows the high energy spec- 
trum of S{q, uj) for U/t = 7.8, V/t = and U/t = 7.8, 
V/t = 3 on a logarithmic scale {L = 60 in both cases). For 
energies higher than the one-particle gap we can see that there 
are dispersive structures up to an energy of roughly twice the 
one-particle gap. The presence of these spectral features fur- 
ther shows that the coupling to the charge sector affects the 
spin dynamics of the system. 



IV. DYNAMICAL CHARGE STRUCTURE FACTOR 

In this section we present numerical results for the dy- 
namical charge structure factor N{q, uj) of the extended Hub- 
bard model. This correlation function is probed by electron- 
electron energy loss spectroscopy (EELS). It is also believed 
to describe the dispersion of peaks and onsets in RIXS exper- 
iments'*. 

In Fig.|5]we show a density plot of N{q, Lo)/q^ for the pa- 
rameters U/t ^ 7.8, V/t = 1.3 relevant for SrCu02. The 
present DDMRG data are significantly better resolved than 
in our previous worlii^. Since limg^o / doj iV(q, w) cx q^ 
we have plotted iV((7, u})/q^ on a logarithmic scale in Fig.|5] 
The momentum resolution is Aq/ir = 1/61 and the energy 
resolution is rj/t = 0.1 for a system of L = 60 sites. The 
most prominent feature is a resonance that disperses down- 
ward monotonically starting at the Brillouin zone boundary 
q = TT. In the strong-coupling limit£ the exciton lies below the 
continuum for all momenta. Here, in contrast, there is also a 
substantial spectral weight below the resonance, which clearly 
lies above the low-energy onset of the continuum shown by 
triangles for q > 7r/2 in Fig.|5] As already found in Ref.fla 
this low-energy onset of the continuum follows the spinon dis- 
persion (shifted by the value of the charge gap), which is in- 
dicated by a line for q > 7r/2 in Fig.|5] 

The presence of spectral weight below the resonance is par- 
ticularity blatant at the Brillouin zone boundary. This can be 
seen clearly in Fig. |6l where we show spectra N{q « tt, uj) 
calculated for 60 and 120 lattice sites and the same broadening 
rj/t = 0.1. Both spectra are almost indistinguishable which 
excludes the possibility that the observed spectral weight is a 
finite-size effect. In the inset of Fig.|6lthe same data are shown 
on a logarithmic scale. We note that more spectral weight lies 
below the resonance than above it. Just above the resonance 
there is a small dip in N{q « tt, uj) reminiscent of an asym- 
metric Fano line shapesifi^ 

In our previous worki^ we reported that the extended Hub- 
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FIG. 6: Charge structure factor at g = tt for U /t — 7.8 and V/t = 
1.3 calculated with a broadening r]/t = 0.1. System lengths are 
L = 60 (full line) and L = 120 (open circles). The inset shows the 
same data on a logarithmic scale. 



FIG. 7: Scaling analysis of the logarithm of the peak height h in 
N{q ~ TTjU}) as a function of the logarithm of the size-dependent 
broadening ri{L) = 12t/L for various interaction strengths. Lines 
are linear fits. 



bard model with U/t = 7.8, V/t = 1.3 and t = 0.435 
eV gives an accurate description of the low-energy optical 
conductivity cr(w) = u!limg-,oN{q,uj)/q'^ of SrCu02. Fur- 
thermore, we showed that the low-enery onset and peak dis- 
persions in N{q,Lu) are in quantitative agreement with simi- 
lar dispersive structures observed in a RIXS experiment. We 
can also compare our DDMRG results for N{q, cu) in the ex- 
tended Hubbard model to the EELS spectrum of Sr2Cu03i5S. 
As the model parameters for this material are close to those 
for SrCu02 (see below), we expect to observe the same qual- 
itative features in both cases. Actually, we find that our re- 
sults for the charge structure factor are in qualitative agree- 
ment with the EELS spectrum of Sr2Cu03. In particular, the 
low-energy tail of N{q, uj) explains the signal observed exper- 
imentally below the resonance close to the zone boundary }2i a 
feature which could not be explained by the strong-coupling 
theory. 

To understand the nature of the resonance observed in 
N{q « TT,u!) we have performed a finite-size analysis of the 
peak height h{ri{L)) for different interaction strengths. The 
broadening is size-dependent such that 7]L = 12tM'ln strong- 
coupling theorySthe dispersive resonance is a (5-peak indicat- 
ing the presence of an exciton. As discussed in Ref. [3 the 
peak height should diverge as h{ri) oc ?7~^ in that case. Fig- 
ure shows a plot of h{ri) on a double logarithnic scale for 
various couplings. Clearly, h{i]) diverges as power law 77" 
for rj ^ 0. We extract the exponents a using a linear fit 
of the logarithms of h{ri) and ?]. The error in the value of 
the exponent a is 0.01 or smaller For U/t < 10 we clearly 
find that > a > — 1, which indicates a power-law diver- 
gence in the spectrum N{q,uj) but no (5-peak,!?. In fact, the 
prediction of the strong-coupling expansion, a « —1, is only 
reached (within the numerical accuracy) for extremely large 
U/t > 100 (not shown in Fig.0. This is a first indication that 
strong-coupling theory becomes accurate only for unphysi- 
cally large on-site Coulomb interaction. 



Since the onset of N{q, uj) follows the spinon dispersion, it 
is clear that the excited states contributing to the charge struc- 
ture factor consist in charge excitations hybridized with spin 
excitations. This is the reciprocal effect of the itinerancy ef- 
fect in the spin structure factor. The Fano-type line shape of 
the resonance suggests that it originates from a discrete charge 
excitation coupled to a continuum of spin excitations. In the 
strong-coupling limit this discrete charge excitation is an ex- 
citon (a neutral bound state of a holon and an anti-holon). As 
the spin excitation band width goes to zero in that limit, only 
elastic scattering is possible and thus the exciton has an infi- 
nite lifetime and appears as a (5-peak in the spectrum N{q, lo). 
For a finite interaction strength, however, spin excitations have 
a finite band width. Thus neutral bound states have a decay 
channel into independent charge excitations due to the (inelas- 
tic) scattering by spin excitations and therefore their lifetime 
is finite. Accordingly, the resonance is no longer a (5-peak but 
a power-law divergence. In particular, for the couplings U 
and V relevant for real materials, such as the Mott-insulating 
chain cuprates SrCu02 and Sr2Cu03, the resonance is not an 
exciton but corresponds to neutral bound states with a finite 
lifetime. 

As a further check of the strong-coupling theory we directly 
compare its predictions with our DDMRG results in Fig. |S] 
The DDMRG data were obtained for a 30-site lattice and a 
quite large broadening r\/t ^ 0.4, which smears out fine de- 
tails. We have convolved the strong-coupling structure factor 
with a Lorentzian of the same width t] to enable a direct com- 
parison. 

In Fig.iJa) we show N{q, for U/t = 7.64, V/t = 2.36. 
These parameters were chosen in Ref. to explain EELS 
data of Sr2Cu03 using the strong-coupling approach. How- 
ever, the comparison in Fig. |8ja) reveals large deviations 
between DDMRG and strong-coupling results for these pa- 
rameters. Therefore, our results demonstrate that a strong- 
coupling approach is not accurate in this parameter regime 
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FIG. 8: Comparison between the strong-coupling- (solid lines) and 
DDMRG (dashed lines) results for the charge structure factor for two 
sets of model parameters. A broadening rj/t = 0.4 is used for both 
approaches. The momentum varies from g = O.OStt (top) to q = 
0.977r (bottom). 



and the parameters U/t = 7.64, V/t = 2.36 determined us- 
ing this approach are not reliable. For comparison, a fit of the 
experimental optical conductivity and exchange coupling to 
DDMRG results for the extended Hubbard model (as done 
in Ref. IT5 for SrCu02) yields the parameters U/t — 7.4, 
V/t = 1.8, and t = 0.42 eV for SraCuOs. We note that 
the apparently small difference between strong-coupling and 
DDMRG parameters has actually radical effects, as the opti- 
cal absorption spectrum contains an exciton below the Mott 
gap for V > 2t but not for V < 2t. (Unfortunately, due to the 
finite experimental resolution it is not clear whether there is 
such an exciton in the linear optical absorption of Sr2Cu03.) 

As expected the agreement between strong-coupling the- 
ory and DDMRG becomes better for stronger interactions. 
This can be seen in Fig. ISjb), where we compare N{q,uj) 
for U/t = 20 and V/t ~ 1. Nevertheless, in spite of the 
interaction equaling five times the bare band width At, we 
still observe significant deviations from the strong-coupling 
predictions. Similar discrepancies were found previously at 
coupling as large as [/ = AOt in a study of the optical con- 
ductivity^^ which is related to N{q,uj) through (t{uj)/u! = 
\iing^oN{q,uj)/q'^. 

It has been shown in Ref. that there are four different 
types of optical excitations that contribute to (t{lu) depending 
on the strength of the next-neighbor repulsion V. We now 
study the evolution of N{q = tt, lu) with V and discuss how 
the nature of these excitations influences the dynamical charge 
structure factor. Figure |9] shows N{Tr,uj) for U/t ~ 7.8 and 
several values of V/t. As expected, we find a broad peak when 
there is no next-neighbor interaction since the charge excita- 
tions do not form long-lived bound states. For V/t = 1.3 we 
find a sharp peak with a small but finite intrinsic width. Wlien 
V/t = 2 the resonance is weaker and spectral weight leaks 
to lower energies which leads to a pronounced asymmetry of 
the peak. This trend continues until the peak is no longer vis- 




FIG. 9: Charge structure factor N{q ~ vr,tj) for U/t = 7.8 and 
V/t = (full line), V/t = 1.3 (dotted line), V/t = 2 (short dashed 
line), V/t = 2.75 (long dashed line), and V/t = 3.5 (dashed dotted 
line). A broadening r) = Q.lt is used. 



ible in the continuum. The spectral weight is now distributed 
in a broad band as seen for V/t = 2.75 in Fig. O In the 
optical conductivity CDW droplets are known to give rise to 
such a bandiS Therefore we suggest that the band observed 
at (7 = TT is also related to CDW droplets. When V is fur- 
ther increased a peak appears at low energies which contains 
most spectral weight. Approaching the CDW phase boundary 
{Vc/t w 3.5 for U/t = 7.8) the total spectral weight in that 
peak iV(7r) = / da;A^(7r,a;) (i.e, the static charge structure 
factor at g = tt) becomes extremely large, which signals the 
occurrence of the long-range CDW order for V > Vc- 



V. ONE-PARTICLE SPECTRAL FUNCTION 

We now consider the one-particle spectral function A{q, uj) 
of the Hamiltonian Q. This dynamical correlation func- 
tion corresponds to the spectrum measured in ARPES exper- 
iments. A density plot of the spectral function is shown in 
Fig-Elfor the parameters appropriate for SrCu02, U/t ~ 7.8 
and V/t ~ 1.3. These data have been calculated on a 90- 
site lattice using a broadening 77 = O.lt, which yields a mo- 
mentum resolution Aq = 7r/91 and an energy resolution of 
the order of 77. With this resolution spin-charge separation is 
clearly seen in Fig. [TO] The spinon branch is denoted by (S) 
and has a width Ws = ■nJ/2 w QM w 0.4 eV. The holon 
branch is labeled (H) and has a width of about 2.7t ~ 1.2 
eV (from q — to q = 7r/2). This spin-charge separation 
has been observed experimentally in the ARPES spectrum of 
SrCu02»2i^ The spinon and holon band widths calculated here 
for U/t = 7.8 and V/t = 1.3 are compatible with the ex- 
perimental data: The latest measurements? yield 0.49 ± 0.13 
eV and 1.42 ± 0.08 eV, respectively. [Note that the effective 
hopping term t defined in Ref. is not equivalent to the bare 
hopping term t of our model Q.] A further verification of 
the theoretical description of SrCu02 by the extended Hub- 
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FIG. 10: Density plot of the one-particle spectral function A{q,u}) 
calculated with DDMRG for U/t = 7.8 and V/t = 1.3. The dis- 
persions of elementary excitations in the one-dimensional Hubbard 
model with U/t = 7.74 are also shown: Spinon branch (S), holon 
branches (H,H',H") and lower boundary of the holon-spinon contin- 
uum (LB). Open white circles mark the position of a step-like feature 
in the DDMRG spectrum. 



bard model would be a direct comparison of the experimental 
spectral weight distribution with our DDMRG data. Unfor- 
tunately, the limited resolution of ARPES data and a strong 
background signal prevent a meaningful comparison. 

In addition to the holon and spinon branches one sees sev- 
eral dispersive structures in the spectral function of the ex- 
tended Hubbard model, which are quite similar to those found 
in the Hubbard model {V ~ 0)iiii^ To identify these struc- 
tures we have calculated the exact dispersions of various ex- 
citations in a periodic 90-site Hubbard chain using the Bethe 
Ansatz. The parameter U/t = 7.74 has been chosen to repro- 
duce the one-particle gap of the extended Hubbard model with 
U/t = 7.8 and V/t = 1.3. The dispersions of the Hubbard 
model corresponding to the most important features in the 
spectral function of the extended Hubbard model are shown 
as Unes in Fig.^| There are the spinon branch (S), the holon 
branch (H), a secondary holon branch (H'), the continuation 
of the holon branch (H"), and the lower boundary (LB) of the 
continuum of excitations consisting of a single spinon and a 
single holon . In the entire Brillouin zone these dispersions of 
the Hubbard model follow closely the dispersive structures of 
the extended Hubbard model. The deviations can be mainly 
attributed to a change of the spinon and holon band widths. 
Thus we conclude that in the Mott-insulating regime the spec- 
tral functions of the extended Hubbard model are qualitatively 
similar to those of the Hubbard model. 

There is a very weak dispersive structure in A{q, oj) of the 
extended Hubbard model which is not visible in the density 
plot and is indicated by open circles in Fig.^| This feature 
is not present in the spectral function of the Hubbard model 
with U/t = 7.74. In the extended Hubbard model it consists 
of a shoulder in the spectrum and can be localized by cal- 
culating the first derivative of A{q, cu). Its dispersion closely 
resembles the holon dispersion shifted by an offset equal to 
the spinon band width. However, we do not find this shoulder 
if we use periodic boundary conditions. Therefore, this fea- 
ture is a boundary effect due to the abrupt cutting-off of the 
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FIG. 11: One-particle spectral function A{q,Lu) for various mo- 
menta q and three different interaction strengths: U/t = 7.8, 
V/t = 1.3 (dotted hne), U/t = 7.8, V/t = 3 (dashed line) and 
U/t = 7.74, V/t = (full line). The broadening is rj/t = 0.1. 



nearest-neighbor interaction terms at both chain ends. 

We note that the spectral function A{q,u!) of the extended 
Hubbard model has been investigated previously using a vari- 
ational CPT approach'^ for a set of parameters {U/t = 
8, V/t = 2) close the one used in Fig.^| Qualitatively, both 
spectral functions are very similar and resemble the spectral 
function of the Hubbard modeIiiiii2ii However, our DDMRG 
spectra are free of the striped structures which are present in 
the variational CPT spectral functions. Therefore, we con- 
clude that these striped structures are artifact of the variational 
CPT approach. 

We now discuss the impact of varying the model param- 
eters on the spectral function A{q,uj). A comparison of the 
line shapes is shown in Fig. ^2 for various momenta q and 
three different parameter sets (U/t ~ 7.8, V/t ~ 1.3), 
{U/t = 7.8, V/t ^ 3) and (U/t ^ 7.74, V/t ^ 0). These data 
have been calculated for a broadening ij/t = 0.1 on a 90-site 
chain. For \q\ < qp ~ Tr/2 there are only small differences be- 
tween the line shapes for V/t ^ and V/t = 1.3. The holon 
and spinon peaks shift according to the change of the holon 
and spinon band widths. This is clearly visible in the two up- 
per panels of Fig.[^ where spinon and holon peaks are well 
separated. For a stronger coupling V/t = 3a significant trans- 
fer of spectral weight occurs from the spinon peak to the holon 
peak. The low-energy continuation of the holon branch (H") 
is also visible as a small third peak in the right upper panel. 
The two middle panels of Fig. ^2 shows the spectral function 
as q approaches tt/2 from below. There, spinon and holon 
structures appear to be merged in a single strong peak because 
the energy difference between spinon and holon is small com- 
pared to the broadening i]/t = 0.1 used. In these panels one 
again sees the peak associated with the continuation of the 
holon branch (H") at high binding energies {us/t w —6). The 
spectral weight of this structure diminishes as V increases and 
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this structure is no longer visible for V/t>'i. 

For momenta \q\ > Tr/2 the total spectral weight is much 
smaller than for \q\ < tt/2, as expected, and diminishes for 
increasing V. Nevertheless we observe several spectral fea- 
tures in the two lowest panels of Fig. ^2 which are usually 
referred to as shadow bandsi^S Shifts of the peak position for 
varying U and V are again related to the change of the spinon- 
and holon band widths. The distribution of spectral weight is 
more strongly affected by the value of the nearest-neighbor in- 
teraction for \q\ > tt/2 than for small momenta q. The peaks 
associated with the shadow band and the secondary holon are 
well separated in the left panel {q < tt) but have merged in 
the right panel (q tt). Furthermore, a low energy struc- 
ture corresponding to the lower boundary of the holon-spinon 
continuum is also visible around cj/t = — 7. 

Finally, we have found that there is a little spectral weight 
below the lower boundary of the holon-spinon continuum. 
This is most clearly seen in the lowest left panel of Fig. ^2 
for V/t = 3. Therefore, excitations made of one holon and 
one spinon are not sufficient to explain all features of the 
one-particle spectral function of the extended Hubbard model. 
Overall, the influence of the nearest-neighbor interaction V 
on the spectral function is noticeable but not as dramatic as 
its impact on the spin and charge structure factors. The main 
features (spinon, holon and shadow bands) are always present 
with similar dispersions albeit with different spectral weights 
and band widths as already found with the variational CPT 
methodii^. 



VI. SUMMARY 

We have calculated the dynamical spin and charge struc- 
ture factors and one-particle spectral function of the one- 
dimensional extended Hubbard model with on-site (U) 
and nearest-neighbor {V) interactions using the dynamical 
density-matrix renormalization group method. We have in- 
vestigated how these dynamical correlation functions are af- 
fected by variations of the model parameters and discussed the 
spin and charge dynamics of this model in the Mott insulat- 
ing regime. For the spin structure factor the strong-coupling 
picture of the Heisenberg model remains qualitatively valid 



down to intermediate on-site repulsion U At, in particu- 
lar for the parameters relevant for cuprate compounds. Itiner- 
ancy effects due to the coupling with charge fluctuations be- 
come significant for weak interactions U < 3t. The nearest- 
neighbor interaction V does not affect the spin structure fac- 
tor qualitatively. For the charge structure factor, however, the 
strong-coupling approach fails but for extremely strong in- 
teraction U ^ t because the coupling to spin fluctuations 
becomes relevant as soon as t/U is finite. In particular, it 
leads to incorrect predictions for the charge dynamics for the 
parameters appropriate for cuprate compounds. Varying the 
nearest-neighbor interaction V causes dramatic alterations in 
the charge structure factor, which reflect changes in the nature 
of the charge excitations, such as the formation of (quasi-)- 
excitons and CDW droplets. The one-particle spectral func- 
tion does not change qualitatively with varying parameters 
within the Mott insulating regime. 

Finally, we have shown that the extended Hubbard model 
allows for a quantitative description of various experiments 
in the cuprate chains SrCu02 and Sr2Cu03 with only a sin- 
gle choice of parameters for each material. For SrCu02 the 
extended Hubbard model with U/t = 7.8, V/t = 1.3 and 
t = 0.435 eV describes the low-energy optical conductivity, 
the main dispersive features in the RIXS spectrum, and the in- 
elastic neutron scattering spectrum quantitatively. In addition, 
the one-particle spectral function agrees at least qualitatively 
with the experimental ARPES spectrum. For U/t — 7.4, 
V/t ^ 1.8, and t = 0.42 eV the extended Hubbard model 
fits the low-energy optical conductivity spectrum and the ex- 
change coupling of Sr2Cu03 very well. Moreover, it repro- 
duces the main features of the EELS spectrum of that mate- 
rial at least qualitatively. In conclusion, the extended Hubbard 
model yields a coherent quantitative description of low-energy 
electronic excitations in cuprate compounds and thus provides 
a theoretical framework for studying the spin and charge dy- 
namics of these systems. 
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